Entanglement of Gaussian states using beam splitter 
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We study an experimental scheme to generate Gaussian two-mode entangled states via beam 
splitter. Specifically, we consider a nonclassical Gaussian state (squeezed state) and a thermal 
state as two input modes, and evaluate the degree of entanglement at the output. Experimental 
conditions to generate entangled outputs are completely identified and the critical thermal noise to 
destroy entanglement is analytically obtained. By doing so, we discuss the possibility to link the 
resistance to noise in entanglement generation with the degree of single-mode nonclassicality. 



I. INTRODUCTION 

In quantum information processing more and more at- 
tention is directed to the continuous variable (CV) sys- 
tems as they have emerged as an alternative resource to 
the discrete level systems. The CV states of considerable 
importance are the Gaussian ones. The interest in this 
special class of states stems from the experimental fea- 
sibility to produce them from reliable sources 01 and to 
control them using accessible tools such as beam splitters, 
phase shifter, and squeezers 0. The two-mode entangled 
Gaussian states have been utilized in many of quantum 
information applications [3, 0, S S 01 • 

Entanglement between two Gaussian modes is rou- 
tinely generated in laboratory, e.g., the two output beams 
of a nonlinear optical device (parametric down converter) 
Alternatively, a beam splitter, one of the linear opti- 
cal devices, can also be used to generate quantum entan- 
glement between two modes Q. There have been many 
studies for producing entanglement using beam splitter 
as an ent angler 

H, [13, [m, [Hill, [II- In particular, Kim 
et al. studied the entangler properties with many differ- 
ent input states, such as Fock states, pure- and mixed- 
Gaussian states. They conjectured that in order to ob- 
tain an entangled two-mode state out of a beam split- 
ter, it is necessary to have a nonclassical state at one 
input, which was later proved in [T5| . Furthermore, the 
sufficiency of single-mode nonclassicality to generate en- 
tangled states via beam splitter was demonstrated by As- 
hoth et al. jl6j . On another side, the separability criteria 
to detect such entangled outputs via beam splitter have 
also been derived [ij, [3 ■ Remarkably, a broad class of 
uncertainty inequalities was presented to detect entangle- 
ment produced using generalized single-mode nonclassi- 
cal states, which include higher-order amplitude squeez- 
ing and high-order photon statistics [isj . 

In this paper, we investigate in detail the entangle- 
ment generated via beam splitter using two uncorrelated 
Gaussian input modes. In particular, we consider a non- 
classical single-mode state (squeezed state) and a ther- 
mal state at two input modes, respectively. We note that 
Wolf et al. also considered a closely- related problem, i.e., 
they derived the conditions to generate bipartite Gaus- 



sian entanglement using passive transformations, focus- 
ing on the optimal scheme [l^. They proved that a 50:50 
beam splitter is the optimal choice regardless of exper- 
imental parameters, and remarkably, that the optimal 
degree of entanglement depends on the smallest eigen- 
values of the input covariance matrix. In other words, 
the degree of entanglement is solely determined by the 
degree of nonclassicality regardless of the purity of input 
state. In realistic situations, however, there always oc- 
curs an experimental error in designing beam-splitter, so 
a careful analysis of nonoptimal cases is further required. 

In this respect, we first want to identify the whole ex- 
perimental conditions to successfully generate entangled 
Gaussian states. Second, a deeper issue is to establish 
the link between single-mode nonclassicality and two- 
mode entanglement in general. A specific question we 
address is whether there exists a monotonic relation be- 
tween the degree of input nonclassicality and the critical 
temperature (degree of noise) at which the output entan- 
glement disappears. For this purpose, we parametrize an 
arbitrary single- mode Gaussian state in terms of the non- 
classicality depth T [13] and purity u, and study the be- 
havior of entanglement as a function of t, u,n ( thermal 
photon number at input), and the beam splitter trans- 
mittance. We demonstrate that the monotonic relation 
between the single-mode nonclassicality r and the criti- 
cal thermal noise rTc exists only at the optimal choice of 
50:50 beam splitter and that Uc is generally a function of 
u as well as r for other choices of beam splitter. 

This paper is organized as follows. In Sec. II, a single- 
mode Gaussian state is briefly introduced with its covari- 
ance matrix in terms of nonclassical depth r and purity u. 
In particular, the covariance matrix of the two-mode out- 
put state via a beam-splitter is obtained for the case that 
a nonclassical (squeezed) Gaussian state and a thermal 
state are used as two input modes. In Sec. Ill, the degree 
of entanglement at the output is evaluated in terms of the 
logarithmic negativity and the experimental conditions 
to successfully generate entangled output are derived to- 
gether with optimal setting. The critical thermal noise 
to destroy entanglement is analytically obtained and dis- 
cussed in relation to the degree of single-mode nonclassi- 
cality. In Sec. IV, our main results are summarized with 
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concluding remarks. 



II. TWO-MODE STATES OUT OF BEAM 
SPLITTER 

Consider a lossless beam splitter whose input ports are 
fed by two single mode fields with complex amplitude 
Qfi and a2, respectively. The complex amplitudes of the 
fields at the output ports are given by 



Pi 
P2 



"1 
a2 



(1) 



where Mb is the beam splitter transformation matrix 
given as 



Mr = 



C0S( 

sin 9t 



sin6'e''^ 

-^v cose 



(2) 



The transmittance of the beam splitter is represented by 
cos^ 9 and the phase difference between the reflected and 
the transmitted fields by Lp. 

Nonclassical Gaussian state — Let the first input 
mode, ai, to the beam splitter be a single mode Gaussian 
state defined by a characteristic function of the form 



X (x) = exp 



Ixt^.x 



(3) 



where x''^ = (a*, ai), and Vi is the covariance matrix 

(a: real, h — |&| e*"^: complex). In Eq. we ignore 
the term linearly dependent on x which describes the 
displacement in phase space. This is because the local 
displacement at each input emerges as another form of 
local displacements at the output two modes so that it 
does not affect entanglement property at all. 

A Gaussian state may be classical (coherent and ther- 
mal states) or nonclassical (squeezed states). A number 
of measures have been proposed to quantifythe degree of 
nonclassicality for a single-mode state pol [2ll . [2^ . [23l . [2^ . 
One of them, which will be used in th is p aper, is based 
on the Glauber- Sudarshan P-function [2^ defined as 



1 



d?ax{ct)e' 



lQ|^-aC*-l-a*C 



(5) 



where x("^) = Tr{_D(Q:)/3} is the characteristic function. 
The P-function renders it possible to express the expec- 
tation values of normally ordered operator functions in 
close correspondence to the calculation of mean values in 
a classical stochastic theory. A quantum state is said to 
have a classical analog if its P-function has the properties 
of a classical probability density. In general, however, the 
P-function may fail to be a probability distribution. A 



quantum state is called nonclassical if its P-function is 
either singular or not positive-definite. 

The integral in Eq. ([5|) may not be evaluated for a non- 
classical state in general. However, a smooth and positive 
definite function that becomes acceptable as a classical 
probability distribution is introduced by the convolution 
transformation of the P-function [201 as 

i?(T,?/) = — /d\e-^l«-"lV(C). (6) 

TTT J 

For a given P-function, there exists a certain value of 
Tm such that the P-function becomes positive-definite for 
T >Tm- The threshold Tm generally takes a value in [0, 1] 
and is regarded as a measure of nonclassicality [l^]. In 
case of a Gaussian state with covariance matrix Vi, the 
condition for the positive definiteness of R (r, rj) becomes 



r-i|/>0, 



(7) 



and thus takes a value in [O, i] . From now on, r 
is used instead of to denote the nonclassical depth, 
and for the covariance matrix Vi in Eq. Q, the degree 
of nonclassicality is given by 



T — max 0, — a + \h\ + - 



(8) 



On the other hand, the degree of mixedness in a pre- 
pared quantum state p can be characterized by its purity 
M = tr (p^) ranging from (completely mixed state) to 
1 (pure states). For a Gaussian state with covariance 
matrix Vi, the purity becomes (26j 



1 



2^/detT4 



(9) 



In terms of the degree of nonclassicality t and the purity 
w, one can thus express the elements of the covariance 
matrix Vi of a nonclassical Gaussian state as 

^ = 4^i^(ib7) + i(^-^^)' 



Although the parameter b is complex, its phase does 
not play any role in the output entanglement [Eqs. (|20|) 
and (j2ip ]. Therefore, only two real parameters, r and m, 
are sufficient to describe a general Gaussian state for our 
purpose. 

Thermal state input — Let the second input mode, 
Q2, to the beam splitter be a thermal state defined as 



Pth = 



(1 



^ I") ("I 



where n is the average photon number 



exp 



(12) 



(13) 
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where fcs is the Boltzman constant and T is the absolute 
temperature. The thermal state is a classical Gaussian 
state with the covariance matrix V2 given by 



v. 



n 
f 1). 







(14) 



with purity Mth = 1/ (2n 

Two-mode output — Having one-mode Gaussian 
state of nonclassicality r and purity u at one port and 
thermal state at the other port of the lossless beam split- 
ter, the characteristic function of the two-mode input 
field can be written as x(q^1!'^2) = exp (— ^w^V^nw) , 
where = (a|, ai, 0:2, represents the complex am- 
plitudes of input modes and Mn = Vi ® V2 the covariance 
matrix [23] ■ On the other hand, the beam splitter action 
[Eqs. (HI) and ^] yields the covariance matrix of the out- 



put characteristic function x iPi, P2) 
[vt = (/3*,/3i,/3*,/32)] as 



exp (-jvtFoutv) 



A C 
Ct B 



(15) 



Here, the 2x2 matrices A, B, and C are given by 

,c2 



A 



a cos f + [71 

b* cos2 9 



i)sin2( 



b cos^ 6 
acos'^6 + {n+ i) sin^ 6* 






r._fa sin^ e+(n + ^) cos^ 9 be-^'^f sin^ 9 

^ 6*e2*^ sin^ 9 a sin^ 9 + {n + ^) cos^ 9 ^ 

C = sm0cos0(^V ^^^^21 j. 

(16) 

III. QUANTITATIVE MEASURE OF 
ENTANGLEMENT 

In this section, we study the degree of entanglement of 
the output two-mode states with the covariance matrix 
in Eq. (fT5)) . A state described by a density operator p is 
called separable if it can be written as a convex sum of 
the product states, i.e., 

i 

where ^ ^ 1 and 'Y^^Pi — 1. Otherwise, it is called 
entangled. A number of schemes have been proposed to 
verify quantum entan glem ent between two modes of the 
field [2^, H^, [13, [Sll. |32||. In particular, it was shown 
that PPT criterion is sufficient and necessary for 1 x n- 
modes bipartite Gaussian states. Instead of using these 
criteria, we consider the quantitative measure of entan- 
glement based on the logarithmic negativity [33| . 

The logarithmic negativity is defined as iV = 
log2| Ip^"^! I, where | | = trV A'^ A denotes the trace norm 



71 



N 




FIG. 1: Entanglement (logarithmic negativity) as a function 
of the thermal photon number n and the beam splitter angle 
6 for a pure-state Gaussian input (u — 1) with nonclassical 
depth r. (a) r = 0.2, (b) r = 0.4, and (c) r = 0.45. 



and p is the partially transposed density operator. For 
a general n-mode Gaussian state, the trace norm is de- 
termined by the eigenvalues of — (V^cr)^, the so-called 
symplectic eigenvalues of the real covariance matrix Vr- 
Here, the elements of the symplectic matrix a are given 
by the commutation relations, [Ra, Rp] = iuafi, where 
Ra {a — \, - ■ ■ , 2n) denotes the canonical variables [l^ . 
For the covariance matrix in Eq. (llSp . the characteris- 
tic equation to evaluate the symplectic eigenvalues under 
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partial transposition becomes [33| 

$4 - {Det[A] + Det[S] - 2Det[C]) + DetiKut] = 0, 

(18) 

where A, B, and C are the sub matrices of the matrix 
Vout in Eq. (fT5|) . Let two positive roots of this equation 
be The logarithmic negativity is then given by (ssj 

iV = max{0,-log2(2^_)} + max{0,-log2(2e+)} 
= max{0,-log2(2^_)}, (19) 

because the condition Det[K)ut] — Det[Vin] — ^^"g|^2^ > 
jg always holds so that the larger root emerges as 2^+ > 
1. 

Using Eqs. (fTSj) and (fT6)) . the negativity can be repre- 
sented in terms of parameters t, u, n, and 0, which turns 
out to be 



N = max{0, - ^log^ [s - ^j S'^ - ^^^^^ 



}, (20) 



where 



S=-[{n-T + l)S+-{n + T)S- cos 461] 



with 



S± 



1 



v?{1-2t) 



± (2rT+ 1). 



(21) 



(22) 



Note that the negativity does not depend on the phase 
shift Lp at the beam splitter. In the following, we study 
in detail the degree of entanglement as a function of the 
experimental parameters. 



In Eq. (in 



A. Optimal beam-splitter 

S takes extremal values at 6 

at 



and 



f . We obtain S = S^^ = ^ + (2!i±lL at 61 = 0, and 



at 



As 



S - Sjl = i^ILhii + 1 - 2r^ . - 4 

Sa^iL are both positive, so is S for the whole range of an- 
gles 0. The logarithmic function in Eq. (|^D|) is a monoton- 
ically decreasing function of S and the negativity there- 
fore becomes maximal at the largest value of S. 

For the case of S- > 0, i.e., ^ ■ 2f^i_2r) > (2" + 1), 
in which the thermal photon number is relatively small, 
the maximum value of S occurs at 9 — j (50:50 beam 
splitter). On the other hand, for the case of „2(i"'_2t) 
(27i-|-l), in which the thermal photon number is large, the 
maximum occurs at = 0, which essentially corresponds 
to no beam-splitter action and leads to no entanglement 
at all. Therefore, we conclude that the optimal choice 
of beam splitter is a 50:50 one regardless of all other 
parameters ( r, u, and n.) 




3n 2 




(b) 



FIG. 2: Entanglement (logarithmic negativity) as a function 
of the purity u and the beam splitter angle 9 at a fixed level 
of thermal noise n for the nonclassical depth t = 0.45. (a) 
n = 1 and (b) n = 4. 



B. Case of 50:50 beam splitter 



In this case, the negativity is reduced to 



TV = max{0, -log2\/(2n + 1)(1 - 2r)}. 



(23) 



The degree of entanglement is thus independent of pu- 
rity, u, and depend only on the nonclassicality, r. From 
Eq. P5|) , the critical thermal noise for the vanishing neg- 
ativity, = 0, is obtained as 



(24) 



1 - 2r 



If n > He, the entanglement at the output disappears. 
Note that the critical value Tic is a "monotonic" func- 
tion of nonclassicality r regardless of purity u. There- 
fore, the resistance to noise, Uc, in generating entangled 
output can be understood as equivalent to a measure of 
single-mode nonclassicality, r. At the maximal squeez- 
ing, T ^ ^, the critical value approaches Uc ^ oo, i.e., 
entanglement persists at any level of noise. 
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C. Case of general beam splitter 

For a general beam-splitter angle 0, the critical value 

Uc is obtained by requiring 25" = 1 + ^'^^^i^^ in Eq. (25). 
Unlike the case of 50:50 beam splitter, the critical noise 
ric and the logarithmic negativity N depend on the initial 
purity u as well as the nonclassicality t of the input state. 
This implies that the interpretation of the critical noise 
as a measure of nonclassicality is not valid for a general 
beam-splitter setting. For example, at the choice oi 9 = 
the critical value becomes ric = 0.75 for t = 0.3 and 
u = 1 (pure-state), and Tie ~ 0.36 for r = 0.4 and u — 0.2 
(mixed-state). In this example, the higher nonclassical 
depth leads to the lower critical thermal noise. 

Case of near-optimal BS: Let us denote the beam 
splitter angle by 6* = ^4^, where 5 is a small error. The 
fractional deviation from the optimal transmittance 1/2 
then corresponds to e = | , and the critical thermal noise 
is found to be 
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(a) (b) 

FIG. 3: The critical thermal noise ric as a function of the 
purity u and the beam splitter angle 6 for a fixed nonclassical 
depth T = 0.4. The right plot (b) shows a magnified view 
over a narrow range of u close to 1. 



IV. CONCLUSIONS AND REMARKS 



"-"T^^i^''^ l-.^(l-2r)0 ' ('^^ 

which shows the dependence on the purity as well as 
the nonclassicality of the input state. In case of very 
small squeezing, r ^ 1, the dependence on the pu- 
rity is negligible as Uc « t(1 — 2e^). On the other 
hand, close to maximal squeezing, r ^ obtain 

In the following, we consider in more detail the case of 
general BS angles. 

Pure-state input: Let us consider the case that the 
nonclassical resource at the input is pure, i.e., u — 1. we 
plot the negativity as a function of beam-splitter angle 
9 and the thermal photon number n for a fixed value of 
T (degree of nonclassicality) in Fig. 1. Obviously, the 
entanglement becomes optimal for the choice of 6* = ^ 
and decreases as the angle deviates from j. In this case, 
one can obtain the analytic expression of critical thermal 
noise as Uc — j-r^r ' '^^i^h is remarkably independent of 
angle 9. Therefore, although the degree of entanglement 
varies with the beam-splitter parameter 9, the entangle- 
ment disappears at the same level of noise ric = i-2t 
regardless of 9. 

Mixed-state input: In general, the degree of entan- 
glement grows with increasing nonclassical depth t and 
purity u. In Fig. 2, the logarithmic negativity is plotted 
as a function of the purity u and the beam-splitter angle 
6* at a fixed level of thermal noise. We see that a success- 
ful generation of entangled output occurs in a broader 
range of angles 9 with increasing purity u. In Fig. 3, 
the critical noise Tic is plotted as a function of purity 
u and the beam-splitter angle 9 for a fixed nonclassical 
depth r = 0.4. As the purity u increases, the distribu- 
tion of critical value ric becomes broader with respect to 
the beam-splitter angle, and it eventually becomes flat 
at It = 1, as argued in the previous paragraph. 



In summary, we have investigated in detail the gen- 
eration of entangled Gaussian states via a beam splitter 
using a single-mode squeezed state and a thermal state 
as two inputs. We have identified the condition to suc- 
cessfully produce an entangled state at the output, and 
evaluated the degree of entanglement as a function of ex- 
perimental parameters, i.e., the nonclassical depth r and 
the purity u of nonclassical source, the thermal photon 
number n, and the beam-splitter angle 9 (transmittance). 

We have established the connection between the crit- 
ical thermal noise ric and the nonclassical depth r, and 
showed that the connection gains a strong interpretation 
only at the optimal choice of 50:50 beam splitter. In 
other cases {9 ^ ^), the critical noise is a function of 
the purity u as well as the nonclassical depth t so that a 
higher nonclassicality does not necessarily lead to a more 
robust resistance to thermal noise. 

It is noteworthy that the critical noise in Eq. to 
destroy output entanglement coincides with the amount 
of thermal noise that can be introduced to the input 
Gaussian state to destroy its single-mode nonclassicality 
(squeezing) in a specialized setting: Suppose one starts 
with a vacuum state |0) to produce a mixed squeezed 
state p that has the covariance in Eq. ^ and the non- 
classicality T. In general, p can be expressed in the 
Kraus-sum representation as p = Mi|0)(0|M|, where 
MiM\ = I. Now, if the initial vacuum state is re- 
placed by a thermal state with the photon number nth as 
p' — ^ ■ MiPtiiM- , it is easy to show that the single- mode 
state p' becomes classical at nth = 1^27' "^hich is none 
other than the critical noise in Eq. (|24p . Therefore, the 
two contextually different noises coincide quantitatively 
in the Gaussian regime. Of course, the result nth = T327 
is valid only for t < i, and thus cannot be readily ex- 
tended to non-Gaussian regime, e.g. Fock-states {t — 1). 
Nevertheless, it seems plausible to have such a relation 
even for non-Gaussian states in a different form. 
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We also note that C. T. Lee attempted to connect the 
nonclassical depth r to the thermal photon number re- 
quired to destroy all nonclassical aspects of the state [13] ■ 
However, the link by C. T. Lee is rather formal, and pre- 
cisely speaking;, the parameter t in Eq. ([6]) is a Gaussian 
noise [35| not necessarily arising from a thermal state: 
When an initial state p is exposed to a Gaussian noise as 

p' = y d'(3P{(3)D{f3)pDH(3), (26) 

where the state p is displaced in phase space by the 
amount of /3 with the Gaussian weighting P(/3) = 

■^e^^^, the P-function of the output state p' is none 
other than the convolution in Eq. ([5]), with the identi- 
fication <T = T. (Note that this Gaussian noise is dif- 
ferent from the noise process mentioned in the preced- 
ing paragraph.) On the other hand, if one mixes the 
initial state with a thermal state with at 50:50 beam- 
splitter, then output state possesses the P-function as 
P(r/) = 2R{V2t]), that is, a scaled distribution of the 



R function of Eq. ([6]). Not to mention that the scaled 
function represents a different density operator, it is also 
known that a scaling transform in phase space does not 
generally map a physical state to another physical one 
[36„ .37]. In contrast, our connection of r to critical ther- 
mal noise has a clear physical meaning in a realistic ex- 
perimental scheme. Our scheme is, however, restricted 
to the class of Gaussian states and it is thus desirable 
to extend the current issues to non-Gaussian states in 
future. 
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